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Abstract. Wc introduce a new computable invariant for strong shift equiva- 
lence of shifts of finite type. The invariant is based on an invariant introduced 
by Trow, Boyle, and Marcus, but has the advantage of being readily com- 
putable. 

We summarize briefly a large-scale numerical experiment aimed at deciding 
strong shift equivalence for shifts of finite type given by irreducible 2x2- 
matrices with entry sum less than 25, and give examples illustrating to power 
of the new invariant, i.e., examples where the new invariant can disprove strong 
shift equivalence whereas the other invariants that we use can not. 



1. Introduction. 

The shifts of finite type (STFs) form an important class of symboUc dynamical 
systems which has fundamental applications in mathematics, physics and computer 
science. The classification problem for irreducible SFTs up to conjugacy is generally 
believed to be undecidable as indicated by the examples of Kim and Roush [9] 
demonstrating the difference between Williams' concepts [15] of shift equivalence 
and strong shift equivalence. Indeed, shift equivalence is decidable [8], but it is the 
more elusive strong shift equivalence which encodes this significant problem, and 
one can no longer hope that these properties are one and the same. 

When attempting to prove that two matrices fail to be strongly shift equivalent 
one has access to a large and very diverse family of invariants developed over the last 
decades, see [11] for a summary. Most of these invariants are efficiently computable 
and comparable as they take the form of algebraic numbers, finitely generated 
groups etc. An important invariant was developed by Trow in [14] and Boyle, 
Marcus, Trow in [4]. This invariant takes the form of the class of a certain ideal in 
a certain integral domain (see below in section 2 for details). We shall refer to it 
as the BMT invariant for brevity. The fact that the BMT invariant is not readily 
computable was the starting point and prime motivation behind the present paper. 

With this invariant as basis we introduce here 2 new invariants: The first is 
defined and computable under a slight technical restriction and is proved to be 
equivalent to the BMT invariant. The second invariant is defined uncondition- 
ally, is possibly weaker than the BMT invariant, but has the advantage of being 
computable. 

More precisely, the BMT invariant takes the form of the class (defined in the 
usual way) of a certain ideal of the ring Z[l/A] where A is a certain algebraic integer. 
Under a technical restriction involving the conductor of the order O := Z[A] our 
first new invariant is defined as a certain element of the Picard group of O. This 
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new invariant is shown to be equivalent to the BMT invariant when it is defined. 
Our second invariant is defined unconditionally as a certain element in the class 
group of the algebraic number field Q(A). This second invariant is weaker than the 
BMT invariant in the sense that equality of BMT invariants implies equality of our 
second invariants, but it is computable by standard algorithms in algebraic number 
theory as implemented for instance in the computer algebra package Magma, [3]. 
This leads to an algorithmic approach to testing this necessary condition which, as 
we shall see, is quite efficient in disproving strong shift equivalence where all other 
invariants fail. 

We have combined this contribution with other tools that are already described 
in the literature to perform a complete analysis of the question of strong shift 
equivalence in a limited universe of SFTs, given by irreducible integer valued 2x2 
matrices with an entry sum less than or equal to 25. Building on a project by O. 
Lund Jensen [10] and using standard database tools we have recorded invariants 
for all matrices in this universe (with the purpose of telling isomorphism classes 
SFTs apart) and concrete strong shift equivalences (with the purpose of identifying 
isomorphism classes). 

The net result of these efforts can be summarized as follows. There are 17250 
matrices in the universe described, and hence 148772625 potential questions of the 
type 'are matrices A and B strong shift equivalent?' We can answer 99.9993769 
% of these questions by this approach, viz. 927 such questions are open. We will 
briefly summarize the methods and results of this project in section 3 below, but 
postpone a full description of it, along with a quantitative analysis of the efficiency 
of the various invariants, to [6]. 

2. New invariants. 

Let 5' be an n X 71 matrix with non-negative, integral coefficients. We call S 
irreducible if for every there is fc > such that the (i, j)'th entry of S'^ is 

positive (S^ is defined to be the identity matrix; in other words, the irreducibility 
condition is empty for the diagonal entries of S) . Irreducibility of S corresponds to 
irreducibility of the associated SET in the sense that any pair of legal words u, v 
can be interpolated by a third word, w, to obtain a legal word uwv. 

Under these conditions, one knows, cf. for instance Theorem 4.2.3 (Perron- 
Frobenius Theorem) of [11], that S has a positive eigenvalue A that occurs with mul- 
tiplicity 1 in the characteristic polynomial of S and whose corresponding eigenspace 
is 1-dimensional, and is such that \fi\ < A for any other eigenvalue fi. Further, the 
eigenspace corresponding to A is generated by a positive eigenvector, i.e., a vector 
with positive coordinates. This uniquely determined eigenvalue is referred to as the 
Perron eigenvalue of S. 

If now A is the Perron eigenvalue of S there is a corresponding eigenvector 

V = (ui, . . . ,t;„) 

with coordinates wi, . . . , w„ in the ring Z[A]. As the eigenspace corresponding to A 
is 1-dimensional, the vector v is uniquely determined up to multiplication with a 
non-zero element of the algebraic number field Q(A). 

We can then define the BMT invariant of S as the class I{S) of the ideal: 

Rvi + . . . + RVn 
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of the rmg R ;= Z[l/A]. Here, 'class' has the usual meaning: Ideals C and D oi R 
are called equivalent if there is ^ G Q(A) such that £,C = D. 

The BMT invariant is an invariant because of the following statement that follows 
from Theorem 12.1 in section 12.3 of [11] (see also Theorem 6.1 of [4], as well as 
[14]): Suppose that S and T are matrices with integral, non-negative entries that 
are irreducible in the above sense and have the same Perron eigenvalue. Then if 
the SFTs attached to S and T are strongly shift equivalent, we have T{S) = I{T). 

In Theorem 1 below we introduce 2 new invariants. The first of these is not 
always defined, but when it is defined for both S and T it coincides for S and T 
if and only if T[S) = I(T), and is in this sense equivalent to the BMT invariant. 
This invariant takes values in the Picard group of the ring Z[A]. 

The second invariant is always defined, takes values in the class group of the 
algebraic number field Q(A), and is weaker than the BMT invariant in the sense 
that T{S) = I{T) implies that the second invariants of S and T coincide. 

Now let us begin to define these new invariants and prepare Theorem 1 below. 
We will work with the following slightly more general setup and notation: 

K : an algebraic number field, i.e., a finite extension of Q 

Ok ■ the ring of algebraic integers in K 

G\{Ok) ■ the class group of Ok 

A : an algebraic integer in K of the same degree over Q as K 

In other words, the assumption on A is that K = Q(A). Then the ring 

O := 1[X] 

is of finite index in Oki i-e., is an order of K . As a reference for the general theory 
of orders in algebraic number fields, see for instance Chap. 1, §12 of [12]. 

In particular, attached to the ring O is a certain ideal of O, the so-called con- 
ductor of O. Ideals of O prime to the conductor have unique factorizations into 
products of prime ideals. Attached to O is the Picard group Pic(O) of invcrt- 
ible ideals modulo principal ones; the Picard group coincides with the class group 
C\{Ok) if O = Ok- The class group C\{Ok) is a canonical quotient of Pic(O). 

If C is an ideal of Ok we shall denote by [C] the class of C in CI{Ok)', similarly, 
if C is an invertible ideal of O, the symbol [C] denotes the class of C in Pic(O). 

We consider the ring Z[l/A]. This ring is in fact the localization 0(m) of O with 
respect to the multiplicatively closed system 

M:={1,A,A2,A3,...} . 

The claim follows immediately once we notice that A £ Z[l/A]: For A satisfies a 
polynomial equation: 

A"4-aiA"-i-h...-ha„=0 
where n := [K : Q] and the ai are integers. It follows that: 

A = -ai - . . . - a„ ■ e Z[l/ A] . 

For ideals in any one of the rings we are considering above, we have the usual 
equivalence relation denoted by and defined by: C ~ D if and only if there 
exists £,€ K'^ such that ^ D. 
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We will now for the remainder of this section assume that we are given two ideals 
21 and 05 of the ring 0(j\/) . Since 0(a/) is a localization of O we know by general 
theory, see [16], Chap. IV, §10, p. 223, that every ideal of 0(a/) is extended from 
an ideal of O, in other words, that there are ideals A and B oi O such that: 

21 = 0(M) -A, » = 0(M) • B . 
We fix such ideals A and B. 

If C and D are ideals of O we employ the usual notation (C : D) to denote the 
fractional ideal: 

(C : D) := {e e A- I C C} . 

Theorem 1. ('ij. Retaining the above notation, we have 21 ~ *B ij and only if there 
are elements x G {A : B) and y G {B : A) such that: 

xy = A*^ 

for some non-negative integer k. 

(ii). Suppose that the ideals A, B, and O ■ A are all prime to the conductor of O 
and let {P^jig/ be the set of prime divisors in O of the ideal O ■ A. 
Then 2t ~ S «/ and only if: 

[A] = [B] mod ([P,] lie/) 

in Pic(O). 

(Hi). In any case, if {Qj}jeJ denotes the set of prime divisors in Ok of the ideal 
Ok ■ A then a necessary condition for 21 ~ 05 is that: 

[Ok-A] = [Ok-B] mod([Q,] | j e J) 

in C\{Ok). 

Proof. We first prove the 'only if parts of (i) and (ii) as well as part (Hi) simulta- 
neously. So, suppose that 21 05. Since K is the field of fractions of O there are 
nonzero elements a,(3 e O such that a2t ~ /305, i.e., such that: 

0(A/) • aA = 0(A/) • I3B . 

Now, O is a Noetherian ring so the ideals A and B are finitely generated O- 
modulcs. Write: 

A = ^0-a„ , B = ^0-hr 
aes reT 
with certain elements S A, 6t- G B, and finite index sets S and T. 

For each a £ S we then have aoa G 0(m) ' PB and so there is an element So- G M 
such that: 

SrjOtaa G (3B . 

Similarly, there is for each r G T an element G M such that: 

trPhr G aA . 

Putting: 

s := So- , t := , 

o-GS tGT 

we conclude that: 

saA C j3B , tpB C aA , 
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and so consequently, xB C A and yA C B if we put: 

/3 a 
x:=t - - , y= s - - . 
a fi 

We have then a; £ {A : B). y ^ [B : A), and xy = st G M so that xy is a 
non-negative power of A: 

xy = A*' 

for some k G Z>o. 

If now additionally the hypotheses of (ii) are fulfilled then the ideals A and B 
are invertible ideals of O. We can then write {A : B) ~ AB~^, and since now 

O-xC AB-^ 

we have 

O-x^AB-^ -U 
with a certain ideal U oi O (namely. U = A^^B ■ Ox). Similarly, 

0-y = A^^B ■ V 

with a certain ideal V oi O. Then UV = O ■ xy = O ■ which shows first that U 
and V are both prime to the conductor of O (since O • A is), and then that their 
prime divisors arc all among the prime divisors {Pi}ig/ of A in O. Hence, 

[U] G ([P.] I I G /) 

in Pic(O); furthermore, as O ■ x ^ AB^^ ■ U, we have 

[A] [B] + [U]=Q 

in Pic(O). We have shown the 'only if part of (ii). 
For the proof of (Hi) observe that we clearly have: 

xOkB C Ok a , yOnA C OkB . 

Since all (nonzero) ideals of Ok are invertible, we can repeat the above ar- 
guments, substituting A and B by OkA and OkB, respectively, to obtain the 
conclusion of (Hi). 

Let us then prove the 'if part of (i): Suppose that we have elements x E {A : B) 
and y G {B : A) such that xy = A*^ for some non-negative integer k. 

As X € {A : B) we certainly have xQ3 ~ xO(m) ■ B C 0(^m) ■ A = %. On the other 
hand, since y G (B : A) and since xy = A*^ is a unit in 0(j\/) we have: 

2t = 0(A/) ■ A = 0(M) -xyAC 0(m) ■ xB = x^ . 

Hence, 21 = and 21 ~ <8. 

Finally, we prove the 'if part of (ii): By assumption there are then non- negative 
integers Vi for i G / such that: 

W A^B-1[PP 

(we can choose the Vi to be negative since Pic(O) is finite). 
Let 

o-x = i[pr 
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be the prime factorization of O • A; by assumption, each rrii is nonzero. Choose then 
fc e N such that k ■ nii > Vi for each i, put Ui := k ■ nii — Vi, and: 

these are ideals of O prime to the conductor, and we have: 

(**) = n = n = o • . 

iei iei 
Now, by (*) and the definition of V we have 

A'^BV = 0-y 
for some y E . Since V is an ideal of O we have: 

y e A-^BV C A-^B = {B : A) . 
Also, [AB-'^] = [V\ = in Pic(O) because of (**); hence, 

AB-^U = 0-x 

for some x G . Since U is an ideal of O (as all Ui are > 0), we see that: 

X e AB-^U C AB-^ = {A: B) . 

Now, O ■ xy = UV = O • A'^ by (**); changing x by a unit of O if necessary 
we then have xy = A'". By the already proved 'if part of (i) wc conclude that 
2l~ «B. □ 

Remarks: In the setting of Theorem 1, the ideals A and B, as well as the fractional 
ideals {A : B) and {B : A) are all finitely generated abelian groups of rank [K : Q]. 
If A and B are given explicitly via generators then generators for {A : B) and 
{B : A) can be computed. 

The fractional ideals {A : B) and {B : A) are in particular finitely generated 
modules over the order O, and 0-module generators can be found if the ideals A 
and B are given explicitly. Thus, the question of solvability of a single equation 
xy = with X <E (A : B) and y E {B : A) reduces to the question of solvability in 
the order O of a single quadratic equation 

f{xi,...,x,) = A'^ 

where / is a quadratic form with coefficients in K that can be determined algorith- 
mically when A and B are explicitly known. 

In [7] it was remarked that the methods of that work show that there is an 
algorithm for deciding a question like this, i.e., the question of solvability of a 
quadratic equation in an explicitly given order of an algebraic number field. 

Hence, condition (i) of Theorem 1 would become an algorithmically decidable 
criterion if one could somehow limit the fc's that have to be considered to a finite 
number. In a sense, such a reduction to consideration of only finitely many fc's is 
what is happening under the favorable conditions of (m) of the theorem, the main 
point being the finiteness of Pic(O). 

The question of whether condition (i) is algorithmically decidable in the general 
case where one or more of the ideals A, B, and O ■ X are not prime to the conductor 
of O is a more complicated question that we will return to elsewhere. 
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3. The experiment. 

As a part of his thesis work [10], Jensen estabhshed a database coUecting sys- 
tematically obtainable information concerning the classification of SFTs given by 
irreducible 2 x 2-matrices with entry sum less than 25. 

In this database, strong shift equivalences were established using brute force 
searches computing elementary equivalences based on 2x2, 2x3, 3x2, and 3x3 
matrices, supplemented by an implementation of Baker's methods ([1], [2]). 

To establish lack of strong shift equivalence, the following invariants are used: 

• The Jordan form of (disregarding the null space of ^, if necessary) 

• The Bowcn-Franks type groups Z"/p(A)Z" where p is one of 

x - 1, X + 1, 2a; + 1, 22; - 1, - X - 1, + x - 1, + 2a; + 1, 
+ 1, a;^ - 1, a;^ + a; + 1, x^ - a; + 1, x^ - 2x + 1, 2x^ - x - 1, 
2x^ - 3x + 1, 2x^ + 3x + 1, 2x2 + a; - 1, 4x2 + 4x + 1, 4x2 _ 43, _^ 

4x2 _ ^ 

• The invariant of Theorem 1 

Perhaps surprisingly, as we shall describe in [6], it turns out that a full arsenal 
of polynomials p is very important for the efficiency of this approach. 

Example 1. Our database contains a multitude of examples for which the invariant 
introduced here is the only means of telling a pair of matrices apart, e.g. 



A = 



"14 


2 


B = 


'13 


5" 









1 


3 


1 



Indeed, the Jordan forms and all Bowen-Franks type groups are identical for these 
two matrices, but in the setup of Theorem 1 (Hi) we get that [Oa--A] G {[Qj\ \ j G 
whereas [Ok'B] ^ {[Qj] \ j £ J). As a consequence of Theorem 1 (Hi), the matrices 
are not strong shift equivalent. 
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